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Abstract 

The dynamic response of a coated inclusion is considered in the context of active cloaking. The active 
cloak is achieved for a coated inclusion in the presence of membrane and flexural waves. In this paper, 
we investigate the design of an active cloak for a coated inclusion in three frequency regimes: the very 
low frequency (monopole dominated) range, the intermediate range, and the higher frequency range in 
which scattering resonances occur. In the first of these ranges, we validate previous work, which resulted 
in a simple mass-compensation design for the monopole scatterer, while in the second and third ranges, 
a combination of the use of an appropriate coating and the appropriate choice of the amplitudes of the 
active cloaking sources is necessary. We show that such cloaking can indeed be effective in the region of 
strong scattering resonances. We give closed form analytic expressions for the required amplitudes of the 
active cloaking sources in the three frequency regions and provide asymptotic estimates and numerical 
illustrations. 


1 Introduction 

Gloaking of objects from detection by waves governed by the Helmholtz equation is nowadays regarded 
as an interesting and important problem (see, for example, IB? E])- The cloaking methods, well-developed 
for the Helmholtz operator, have also been extended to advanced structured media in electromagnetism, 
elasticity, as well as waves in fluids (see m [lo]). The cloaking system can be understood as a by-pass 
structure, which routes a probe wave around an object, reducing the scattering by the object to a minimum. 

We make a distinction between passive and active cloaking. Passive cloaks do not change their 
properties and remain unchanged for all frequencies and all types of incident waves. The drawback of a 
passive cloak is that its practical implementation can only be effective over a limited range of frequencies 
(see, for example, m, [12]). Active cloaks on the other hand can easily overcome this limitation in their 
effective frequency range, and can also adapt to various types of incident waves. Active cloaking relies on 
being able to radiate compensating fields which cancel out perturbations of the incident or probe field by 
the object to be cloaked [13]. 

In a set of papers na-ini, a continuous model of an active cloak for the Laplace and Helmholtz 
equations was developed. The model was set as a variational problem and was analysed mathematically 
in depth, and the effectiveness of the active cloaking was demonstrated. Einite cloaking devices in the 
continuous framework were introduced to shield any objects within the given “shielded” region. 

In a recent paper [18] , an analytical model of an approximate cloak for flexural waves was developed by 
introducing a finite number of active sources, whose complex amplitudes were chosen to suppress the scattered 
field resulting from an interaction of an incident wave with a finite obstacle. The multipole representation 
was used for the solution, and the required number of multipoles in the scattered field were constrained to 
be zero. Although the truncation of the multipole representation introduces an approximation error, this 
approximation was shown to be entirely adequate for a sufficiently wide range of frequencies. 

In this paper, we investigate passive scattering reduction and active cloaking in the context of both 
the Helmhotz equation (applying to waves in elastic membranes) and the fourth-order Kirchhoff plate theory 
(applying to flexural waves in thin elastic plates). We study these phenomena in three frequency ranges: in 


I 


the first the frequency is sufficiently low, so that wave scattering is dominated by the monopole term, for 
both membranes and Kirchhoff plates. We show that the following mass-compensation formula 

_ Pe ~ (^c) /-. -| N 

1 - (aiK)2 ’ ^ ’ 

is applicable to reduction of scattering by an appropriately chosen coating around the inclusion. Here 
PiiPciPe are mass densities of the inclusion, coating, and ambient matrix, respectively, and ai^ac are the 
interior and the exterior radii of the inclusion coating. This classical formula m has been recently used [20] 
in the context of scattering reduction for coated inclusions at low frequencies for which the dominance of the 
monopole term guarantees the applicability of the criterion 0 for both the membrane and Kirchhoff plate 
waves. We show here that, while this formula is effective for membrane waves in reducing the scattering 
from a second-order to a fourth-order effect, the same is not in general true for Kirchhoff waves, where for 


only very specific choices of material parameters, does the formula (1.1) deliver the reduction of scattering 
to fourth-order. The failure of 0 to deliver the reduction to fourth order scattering occurs because in 
general, an appropriate balance between the Helmholtz and modified Helmholtz parts of the Kirchhoff wave 
expansion requires suitable choices for flexural rigidities, not just densities. We also show that, because 
scattering in the low frequency range is very weak, active cloaking is not necessary. 


As the frequency increases dipole and higher order terms become significant. We show that the choice 
of coating parameters is not sufficient to generate effective scattering reduction in this intermediate frequency 
range. A small number of control sources then delivers effective cloaking. 


As the frequency increases further, we enter the resonant scattering region. These resonances result 
in strongly enhanced scattering amplitudes, which change rapidly with frequency. This causes difficulties to 
arise for the application of active cloaking techniques, as suggested by the following argument. 


One strategy by which active cloaking might be defeated would be for the probe beam to be swept in 
frequency, so that the compensating fields would also have to be frequency-swept in corresponding fashion. 
A circumstance under which this would be difficult occurs when the object to be cloaked has resonant 
scattering properties. In the presence of resonances, it is difficult to avoid time lags between the change of 
the probe frequency and the corresponding changes required in the compensating fields around the cloaked 
object. These time lags inevitably increase with the sharpness (or Q factor) of the resonance. Thus, for 
active cloaking to work well near resonant frequencies, one needs to be able to modify the configuration 
of the cloaked object, to avoid rapid changes in scattering properties with frequency. The configurational 
change we investigate here is the use of the coating around the inclusion to deliver a flat reflectance as a 
function of frequency, rather than a very low reflectance. (The latter is inevitably narrow-band at any but 
the lowest frequency range.) This we show enables the active cloaking of the coated inclusion to operate 
effectively over the frequency region of flat reflectance. Up to twelve active sources may be required, but 
good quality cloaking is always possible by this method. 

The paper is structured as follows. Section includes the governing equations and the multipole 
representations of flexural waves near a coated inclusion. An illustration of resonance regimes and active 
cloaking by point sources is given in Section for the case of membrane waves, as a simple precursor to 
the more involved discussion necessary for flexural waves. In Section]^ it is shown that active cloaking is 
still effective in resonance regimes at higher frequencies, as a result of the strategy described above of using 
an appropriate coating to move the inclusion resonances into the desired frequency range. A closed-form 
expression is given controlling the choice of inclusion and coating parameters, which delivers a zero monopole 
scattering coefficient. This generalises the classical result of Konenkov [21] to coated inclusions and active 
cloaking. The Supplementary Material gives a full analytic description of our cloaking method for both 
membranes and Kirchhoff plates. 


2 Problem formulation: scattering of flexural waves by a coated 
inclusion in a Kirchhoff plate 

The out-of-plane elastic displacement (x; t) in a Kirchhoff plate satisfies the equation of motion 

= 0, x = {xi,X 2 ) € 0,k, k = i,c,e, (2.1) 
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where represents active point sources placed in the exterior of the coated inclusion; the explicit 

representation for T and positioning of the point sources are discussed in Section |3.4[ We note that the 
super/sub-script k is to be replaced by i,c or e, depending on whether we are considering the inclusion, 
coating or exterior, respectively. In equation (2.1), is the biharmonic operator, a dot on the variable 
denotes the derivative with respect to time t, p is the mass density, h is the thickness of the plate, is the 
corresponding region, Dk = /[12{1 — z^|)] is the flexural rigidity, with Ek the Young’s modulus and Uk 

the Poisson’s ratio of the corresponding elastic material. 

Assuming time-harmonic vibrations, that is (x) exp(icc;t), outside the support of the 

function *F(x;t) the governing equation (13 can be reduced to 

2 pkhLd‘^^ 


Dk 


w 


(k)( 


(x) = (A+ {A- I3l) = 0, X G (2/,, k = i,c,e. 


( 2 . 2 ) 


Here P‘1 = pkh/Dk is the spectral parameter. Consequently, can be written as the sum of 

solutions of the Helmholtz and modified Helmholtz equations 


= wP{x) + w^^\x), 


(2.3) 


where 


(x) + (x) = 0, (x) - (x) = 0. (2.4) 

In the case of a circular geometry, the inclusion, the coating and the exterior are defined as follows: 
Vti = {^ \ x\-\- x\ < a^}, = {x : + ^2 < Cle = \ Qc, 

as shown in Fig. 


(2.5) 
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Figure 1: A plane wave wo{xi) is incident on the coated inclusion along the xi-axis. The inclusion, coating 
and exterior are associated with the subscripts i, c and e, respectively. 


On the interface boundaries r = ai where the inclusion and coating meet, and r = ac where the 
coating and exterior meet, the conditions for perfect bonding, that is the continuity of the displacement, its 
normal derivative, the moment and the transverse force, lead to 


'U;C)(x) = 


dr dr 


( 2 . 6 ) 


and, with spatial dependence suppressed, 

Uc f dw^^^ 1 \ 

c ^ g ^2 ^ ^ r d0‘^ ) 


= Di 


dr‘^ 


ui f dw^^^ 1 \ 

r \ ^ r d0‘^ ) 


(2.7) 


Dr 


d ('c) I — r'c f dw^^^ 
dr ^ d0‘^ \ dr 


w 


(c) 


= Dj 


On r = aj we have I = i and on r = we have I = e. 


d 

dr 




+ 


\-yid^ 

r2 QQ2 


/ dw^^^ 
\ dr 



( 2 . 8 ) 
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2.1 Multipole representation for a plane wave scattered by a coated circular 
inclusion 


Excluding the positions of point sources, the fields w^’'\ and can be written 

CX3 

(r, 0)= 


AnO 


(2.9) 


where k is either i,c or e, depending on the respective region. We note that when k = the coefficients 


En^ and Fn^ are identically zero to ensure that there are no forces or moments at the origin. Using the 




representations (2.9), equations (2.6)-(2.8) give us two matrix relations 




\ 

d(c) 


— ^(ci) 


(0 

n 

(d 


on r = aj 


( 2 . 10 ) 


where and are 4x4 and 4x2 matrices, respectively; and 


aL‘> \ 


^(ec) 




e; 

B. 


(e) 

n 

(e) 

n 

(e) 


_ g(ec) 


eA / 


(c) 

n 

(c) 

n 

(c) 


on r = ar 


( 2 . 11 ) 


V / 


where and are both 4x4 matrices. The representations of these four matrices are given in 

the Supplementary Material. At this stage, it is important to emphasise that reducing the matrices 
and A^^^^ to block-diagonal structure is essential to pursue the derivation of an analytical formula for the 
scattering matrix S. We note that in what follows denote the required matrices in block- 

diagonal form and are the associated matrices, with all these matrices being given in the 

Supplementary Material. 

We use the four matrices and to construct two further matrices 

and as follows: 

.( 2 . 12 ) 




- 1 . 


where the expressions which would represent the results of the matrix inversions in the previous equation are 
too cumbersome to be given explicitly. The final step in the construction is to cascade via a matrix product 
the progression from the interface between the inclusion and the inner coating to the interface between the 
outer coating and the matrix, giving the explicit transfer matrix which incorporates the interface 

conditions at both interfaces: 


^(ei,*) _ ^(ec,*)^(ci,*) 


(2.13) 


The transfer matrix can then be used to give the wave coefficients An \ Bn \ En^ and Fn^ ir 

the exterior region in terms of the wave coefficients An^, Bn^ in the inclusion region: 


/ \ 

D(e) 

-On 

Fi^^ J 





(2.14) 


The final step in this procedure is to eliminate the wave coefficients in the inclusion region to determine 
the 2x2 scattering matrix Sn which relates outgoing wave coefficients in the exterior region to incoming 
wave coefficients in the exterior region. The relevant equations resulting from the procedure are described 
in the Supplementary Material, and lead to the following matrix equation 
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where on the left-hand side we have outgoing wave coefficients, and on the right-hand side incoming wave 
coefficients. If we now assume that the scattering occurs from a plane incident wave travelling along the xi- 
axis, the coefficients Ah\ Bn^ associated with the inclusion, An \ En \ Fn^ a ssocia ted with the coating 


and En\Fn^ associated with the exterior, respectively, can be readily found from (2.14) (see Supplementary 
Material). Note that in this particular case An^ = E 


since 


wo{xi) 


Y, i^Jn{Per)e 


inO 


(2.15) 


and the Bn^ are identically zero. 

It can be shown that the general expression for the scattered field is (see, for example, m, m, m) 

oo 

w(sc){r,0) = Koi/Ser) + y] + 5^^ii:„(/3er)] cos (n6>), (2.16) 

n=l 


where 5*-^, i,j 


p(e) 


coefficient in the right-hand side of the expression (|2.16|) is in fact equal to the scattering coefficient 

l(e) 


= 1,2, denote the entries of the scattering matrix. At this point, we emphasise that the 


Q (using Aq^^ = 1 and Bq^^ =0)- As has been claimed previously 


with a suitably designed coating 

0. (In fact, Eq^^ can be made equal to zero exactly, with an appropriate choice of the material 
parameters of the coating - see (3.5), (3.6) below.) 


2.2 An algebraic system for active source amplitudes in the presence of an 
incident plane wave 

Now, we assume that a finite number of active sources are applied in the exterior of the coated inclusion 
in order to reduce the scattering of an incident plane wave by the inclusion. Here we derive formally an 
algebraic system for the active source amplitudes, which are computed to make the required number of 
multipole coefficients vanish in the representation of the outgoing scattered field. 

In order to find the total wave scattered from the coated inclusion, we identify two model problems 
illustrated in Fig. These problems deal with scattering from the coated circular inclusion, with the centre 
at the origin. In the first model problem, the incident field is a plane wave advancing along the xi-axis, 
whereas in the second problem, the incident field is generated by a point source of unit amplitude placed 
outside the coated inclusion. 


X2 



X2 



Figure 2: Model problems: (a) plane wave incident on a coated inclusion; (b) point source near the inclusion. 


The plane wave problem. The first model problem involves scattering of a plane wave from the 
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coated inclusion. The displacement field satisfies the equation 


=0 in \ P, 


(2.17) 


supplemented by interface conditions ( 2 . 6 )-( 2 . 8 ) on the interior and exterior interfaces of the coated inclusion. 
Here V is the region occupied by the coated inclusion. The total field can be written as the sum of the 


incident field and the scattered field '^[ 5 ^)’ 


(2.18) 


where is the same as defined in (2.15). 


Outside the coated inclusion, the outgoing scattered field has the representation 


(.c) 


(p) _ 

^{sc) - 


00 

[4^) (1) (/3er) + {per) 


AnO 


(2.19) 


For Per ^ 1, the evanescent term in the scattered field can be neglected since Kn{Pe'^) = 0(exp(—/3er)/(/der)) 
(see (9.7.2) in |23]), and hence 




( 2 . 20 ) 


The single source model problem. The second model problem involves scattering of a radial wave 
from the coated inclusion (see Fig. [^b)). The displacement field, solving this model problem, is denoted by 
and satisfies the equation 


_ Q in \ P, 


( 2 . 21 ) 


where (5(x — denotes the Dirac delta function, centred at a.^^\ It is supplemented by the transmission 
conditions (2.6)-(2.8). The solution readmits the representation 


(^’^■)(x)=«;jyyx)+«;|:;t(x)=G(x - a(^'))+«;|:;t(x) 


(sc) 


(sc) 


and hence. 


t^(*.j)(x)~ ^ ji'=l,...,V as r ^ oo, 


( 2 . 22 ) 


(2.23) 


where exponentially small terms satisfying the modified Helmholtz equation, are not shown. Here G(x — 
denotes the Green’s function for the biharmonic operator, is the scattered field due to the unit source 

l^sc; 

at and are constant coefficients. 

Assuming that the coefficients En^ and are given, active cloaking is achieved by introducing a 

set of N control sources of complex amplitudes Bj placed at the points a^-^^ in the exterior of the scatterer 
V. After the truncation to order L in the expansions (2.19), ( 2 . 22 ), we choose A = 21/ + 1, so that the total 
displacement field W(^totai) is approximately equal to the incident field , that is 


N 


= in^E) 




i=i 


w 


is,j) 




(2.24) 


To find Sj, we substitute (2.19) and (2.22) into (2.24), and obtain the following system of linear algebraic 
equations 


E 


(p) 


2L+1 


^ BjA]!^^^ =0, A: = -!/,...,!/. 

1=1 


(2.25) 


This algebraic system holds not only for the flexural waves in Kirchhoff plates, but it is also valid for 


the membrane waves governed by the Helmholtz equation. Of course, in the latter case the coefficients E- 


(p) 
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and will change accordingly. It is assumed that the positions of active sources are chosen so that the 

system (2.25) is non-degenerate. 


In the text below, this algebraic system is implemented for the case of 4, 8 and 12 sources positioned 
on a circular contour of a sufficiently large radius a^, and with its centre coinciding with the centre of the 
coated inclusion. This leads to an approximate active cloaking, which is efficient in the required frequency 
range. The numerical examples, based on our semi-analytic method, are also presented and demonstrate 
clearly the effectiveness of the method even in the resonant regimes. 


3 Resonant scattering of membrane waves by a coated inclusion 

The ultimate aim of this paper is to cloak an inclusion with a suitably designed coating from flexural 
waves in resonant regimes using active control sources. In such regimes, a small change in frequency results 
in abrupt variations in scattering and thus it becomes important to eliminate frequency sensitivity of active 
sources. With this in mind, as a preliminary problem, it is instructive to study the interaction of flexural 
waves in a membrane containing a coated inclusion. This is because solutions to the biharmonic equation are 
superpositions of solutions to the Helmholtz (a mixture of propagating and evanescent waves) and modified 
Helmholtz equations (all of which are evanescent waves). 

In this section, we give an illustrative example for active cloaking in resonant regimes in a membrane, 
where waves are governed by the Helmholtz equation. The out-of-plane displacement of the membrane 
satisfies the following equations and transmission conditions in three regions, with and u^'^\ denoting 

the displacements in the exterior, coating and inclusion respectively: 


(x) + (x) + /(x) =0, X e 


ki 


(3.1) 

where Pk is defined as pk = pk/l^k-, k runs over e, c or i, and the domains Qk are as in (2.5); /i/c, k = e, 

are the shear moduli in respective regions, and /(x) is the force term representing the active sources. 
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Figure 3: A plane wave wo{xi) is incident on the coated inclusion along the xi-axis. The inclusion, coating 
and exterior are associated with the subscripts i, c and e, respectively. 


We denote, as before, 
Uc, respectively (see Fig. [^. 


the radii of the inclusion and the exterior boundary of the coating by and 
Thus, the interface conditions can be written as 


and 






Mi 


Me 


dr dr 

dr dr 


on r = ai^ 


on r = Uc- 


(3.2) 


(3.3) 
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Excluding the positions of point sources, the associated solutions of equations (3.1) are 


vS^\r,e) 


£ [4'=) J„ (/3er) + 4"^ WeV) 


e*”^ 


k = e, 


(3.4) 


(i) 

where En ^ are identically zero. 


3.1 Derivation of the relation between coefficients En'^ and A 


l(e) 


Using the representations (3.4) with the interface conditions (|3.2|) and (|3.3|), it is straightforward to 


as 


derive the outgoing wave coefficient En^ in terms of the coefficient An^ (characterising the incoming wave) 

Ei^^=TnAi^\ (3.5) 

MnMl^ - M^J 


where 


Tn = - 




(3.6) 


Here Tn stands for the scattering coefficient, and the representations for = 1,2) and J\fn are given 

in the Supplementary Material. 


We note that once the relation (3.5) is established, the coefficients An\ En^ associated with the 
coating, as well as the coefficient An^ associated with the inclusion can be readily found (see Supplementary 
Material). 

Also note that the scattered field in this particular case can be represented as follows (compare with 


(2.16)) 


Uisc) = ro4'^(/3er) + ^2i"r„lf4(/3er) cos(n0). 


(3.7) 


n=l 


3.2 Resonant scattering by a high-contrast inclusion with no coating 


In the case of an inclusion with no coating it can be shown that formula (|3.5|) reduces to 


4") = 


J, 






In particular, for the monopole term (n = 0) 

7raf \pi - Pel 


\E, 


(e). 


l^e 


uT as Piai.Pedi ^ 


(3.8) 


(3.9) 


where Piai^ PeCii are non-dimensional quantities (considering a plane wave travelling in the xi direction as 
before, Aq^^ = 1). We also note that the monopole term for static problems corresponds to a net force 
produced by an elastic inclusion, and is consequently zero; this differs significantly from the dynamic case, 
when the inclusion possesses a non-zero mass density. In the dynamic case, the coefficient Eq^^ of the 
monopole term is defined by the formula (3.8), and plays an important role in scattering since it determines 


the scatte ring cross-section (see Chapter 13 of Born and Wolf [24]). For low frequencies, as evident from 
equation (3.9), the coefficient \Eq^^\ is small, but, in general, it is a rapidly varying function of frequency, 
and may indeed show strongly resonant behaviour for high-contrast inclusions (here caused by the ratio of 
Pe/l^i = 10, see Fig. [^. 
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Figure 4: Monopole amplitude for membrane wave scattering by an uncoated inclusion with the parameters 
Pi = 1.5, Pe = 1.0, Pi = 0.1, Pe = cii = 0.50: (a) as a function of cj, (b) Detail from (a) (solid blue 

rectangle) showing (dashed green) and the empirical fit 0.10(solid grey). 


Figure l^a) shows the monopole coefficient \ as a function of uj for a specific choice of high-contrast 
inclusion; whilst the densities of the exterior and the inclusion are comparable {pi = 1.5, pe = 1.0), the shear 
modulus of the exterior is tenfold that of the inclusion {pi = 0.1, /ie = 1-0). A prominent feature of the 
low frequency region (enclosed by a solid blue rectangle), where is small, is that both the function 

itself and its derivative tend to zero with uj (for which a diagram of the scattered field is given in Fig. 
I^b)). Sufficiently away from the low frequency region, the major features of Fig. [^include a strong peak 
near uj = 1.2, a narrow zero near uj = 2.62 and a rapid rise to a maximum near uj = 3.0. The region of 
rapid variation of |F’o^^| with uj between 2.62 and 3.0 would give rise to difficulties for active cloaking of a 
frequency-swept incident wave, and thus will be discussed in detail below. 


The influence of \ on the scattered field is illustrated in Figs. 5(a), 5(b) and 5(c) Fig. 5(a) 


shows 

the total displacement for a small frequency value of uj = 0.3, which clearly shows a low scattering pattern. 
Figs. |5(b)| (with uj = 2.62) and 5(c) (with uj = 3.0) contrast the total displacement field patterns around 
the inclusion for the frequencies corresponding respectively to the first non-trivial zero of and its 

subsequent maximum. In Fig. |5(b)| scattering is primarily due to terms of dipole order, but the wave fronts 
behind the inclusion are practically straight. In Fig. 5(c)[ both monopole and dipole terms are contributing 
to the distortion of the wavefronts behind the inclusion. The scattering patterns change dramatically even 
though the frequency difference is small. 
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(c) 


Figure 5: Membrane wave scattering by an uncoated inclusion for the parameters pi = 1.5, pe = 1-0, Pi = 0.1, 
Pe = 1-0, ai = 0.50. (a) Total field for uj = 0.3, (b) Total field for uj = 2.62 where vanishes, (c) Total 
field for uj = 3.0 where has nearly a maximum. In the last frame of each row, we present a zoomed-in 

version of the displacement inside and in the close vicinity of the inclusion. The first two frames in the 
inset (b) represent the same field outside the inclusion, but with different colour maps, where the first frame 
excludes the internal field, whereas the second frame includes it. The same comment applies to the first two 
frames of the inset (c). In all contour plots horizontal and vertical axes are xi- and X 2 -axes, respectively. 


3.3 Perturbation of a plane wave by a coated inclusion 

We now discuss two ways of scattering control by an elastic coating placed around an inclusion. 
First, we consider the very low frequency regime dominated by the monopole term, leading to the mass- 
compensation criterion for the choice of the parameters of the coating around the inclusion, as mentioned 
in the Introduction (see formula ( |I.I[ )). Second, we analyse the higher frequency regime, which may include 
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resonance modes. 


3.3.1 Reduction of the monopole term for small frequencies 


Formula (3.5) leads to 


\E^ 


(e)| _ 


= ^\pia^ + Pc{al-a^) - peal\LO^ + 0{uj^) as piUi, Petti, pc^c, PeCic ^ 0. (3.10) 

4/ie 


This suggests the design 0) for a coating which compensates the density difference between the inclusion 
and exterior, so that the scattering cross-section at lo w fr equencies is minimised accordingly. In particular 


~ 0 in the low frequency regime implies that (1.1) holds, that is, the average mass density of the 
inclusion combined with the coating equals the mass density of the exterior material. When pc is chosen 
according to the formula 0, whilst the shear moduli of the coating Uc and the exterior pc are chosen to 
be equal to each other, the graph of versus uj is given in Fig. |^a). Comparing this with Fig. j^a) 

and their respective detailed views in Figs|^b) and|^b), it is evident that in the low frequency regime, the 
scattering cross-section is substantially reduced. However, it is also shown in Fig. |^a) (see for example, the 
region enclosed by the red rectangle) that for higher frequencies we ob serve resonance scattering associated 
with rapid variation and increase of Consequently, the formula (1.1) is of no use in such regimes. 




(e)| 






Figure 6: Monopole amplitude for membrane wave scattering by a coated inclusion with parameter values: 
pi = 1.5, pe = 1.0, pc = {pe - Pi(«i/<^c)^)/(l - {di/GicY) ~ 0.635462, Pi = 0.1, Pc = 1.0, Pc = 1.0, ai = 0.50, 
Uc = 0.77: as a function of uj. Note, the curve flattens approximately for uj G [3,4.2]. (b) Detail from 

(a) (green rectangle) showing (orange/dot-dashed) and the empirical fit 0.15cc;^ (grey/solid). 


3.3.2 Passive control of scattering by an appropriate coating 

The method of active cloaking that will be discussed in the next section, has the advantage of being 
efficient even at higher frequency regimes. We illustrate here the limitations of passive design, where the 
parameters of the coating are chosen to control the monopole term in the scattered field. At a frequency 
when the monopole coefficient \ is small but rapidly varying, active cloaking requires a rapid adjustment 
of the source amplitudes (see the red box enclosing uj ^ 2.62 in Fig. [^a)). In these particular regions the 
scattered field increases significantly and thus we use the term “resonant scattering”. To overcome this 
difficulty, the design of the coating can be revised to reduce the variation of the coefficient \Eq^^\ for this 
frequency range. In Fig. Ba), we combine data from Figs|^a), |^a) together with resulting from 

the revised design. The revised choice of parameters gives a “flattened” curve near the resonant scattering 
frequency uj ~ 2.62 (see the region enclosed by the red box). A magnified view is given in Fig. [^b). Figure 
l^c) shows the behaviour of | in the low frequency regime. It is noted that the mass balance criterion for 
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the coating of the inclusion clearly reduces the monopole coefficient from being to 0(cc;^). However, 

this is of little import since the monopole amplitude is already very small even for an uncoated inclusion. 






Figure 7: (a) Monopole amplitude for membrane wave scattering by (i) an uncoated inclusion with parameter 
values: pi = 1.5, pe = 1.0, pi = 0.1, pe = 1-0, = 0.50 (green/dashed), (ii) a coated inclusion with 

parameter values: pi = 1.5, pe = 1.0, pi = 0.1, Pc = 1.0, Pe = 1.0, = 0.50, Uc = 0.77 and pc = 

{pe — pi(ai/ac)^)/(l — {ai/acY) ~ 0.635462 (orange/dot-dashed) and (hi) a coated inclusion with parameter 
values: pi = 1.5, pc = 0.81, pe = 1.0, pi = 0.1, pc = 0.11, pe = 1.0, = 0.5, Uc = 0.77 (blue/solid). 
Note, the red box highlights where the blue/solid curve flattens approximately for uj G [2.4, 2.9], a lower 
frequency range than that seen in the green/dashed and orange/dot-dashed curves, (b) Monopole amplitude 
for membrane wave scattering for the parameters of (i)-(iii) in (a), but for a range of uj highlighted by the 
red box in (a). We see strong variation in for an uncoated inclusion in this frequency range, but it is 

clear that with an appropriate choice of parameter values for the coating, we can create a region where \ 
remains fairly flat, (c) Monopole amplitude for membrane wave scattering for the parameters of (i)-(iii) in 
(a) where uo is small. The upper two curves (green/dashed and blue/solid) in (c) correspond to second order 
variation of \Eq^^\ with uo while the third (orange/dot-dashed) corresponds to a fourth order variation. 
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Figure 8: (a) Total displacement field for membrane wave scattering by a coated inclusion with parameter 
values: pi = 1.5, pe = 1.0, pc = (pe - pi{ai/acY)/{I - {ai/acY) ~ 0.635462, pi = 0.1, pc = 1-0, pe = 1.0, 
ai = 0.50, Uc = 0.77 and cj = 0.3. (b) A zoomed in version of the displacement inside the inclusion, the 
coating and in the close vicinity of the coating. 





Figure 9: Membrane wave scattering: absolute value of the component n, En^ Hn\l3eO's)j of the displacement 
field, where is the radius of circle on which the control sources are placed, as a function of n for cj = 2.62. 
Parameter values as for (blue/solid) curve in Fig. Pi = 1.5, pc = 0.81, pe = 1.0, pi = 0.1, pc = 0.11, 
Pe = 1.0, ai = 0.5, Uc = 0.77, = 1.57. 

Next, in Fig. we present contour plots of the membrane wave amplitude for a small frequency value; 
the inclusion has a coating whose density is chosen according to the formula The parameter values, 

with the exception of the mass density and elastic modulus of the coating, are chosen exactly as in Fig. |5(a)[ 
we choose the elastic modulus of the coating equal to that of the exterior and the frequency value as cj = 0.3. 
As in Fig. [^c), it is apparent that the scattering has been reduced compared with Fig. |5(a)[ however, the 
benefit of such a reduction is very low, as even without the coating there is no significant scattering. 

For the higher frequency near uj = 2.62, we have chosen the coating to “flatten” the monopole term 
curve over the resonant range of oo. We then observe in Fig. [^that the multipole held components of 
higher order become non-negligible, and hence passive reduction of scattering would not be appropriate any 
longer. Active cloaking by point sources, as described in [18], enables us to make the multipole coefficients 
vanish to the required order. 
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3.4 Active sources: cloaking of membrane waves in resonant regimes 

As in m, we add active sources of membrane waves to annul a number of multipole coefficients of 
the propagating scattered field. Our study in [18] confirmed that for an effective cloaking of flexural waves 
in a Kirchhoff plate at least four sources were required. Thus we first attempt to create a cloak for the 
coated inclusion with four active sources. These sources are symmetrically located on the xi- and X 2 -axes, 
all at a distance > Uc away from the origin. Note that the symmetry of the inclusion dictates that we set 
the amplitudes of the sources on the X 2 -axis to be the same. Since the Green’s function for the Helmholtz 
operator is given by 

G{xi -U,X2-v) = (a:i - + {X2 - v)"^), 

where (u^v) is the position of the point source, using Graf’s addition theorem for the sources located at 
(—0), (ttg, 0), (0, —as) and (0, a^) we obtain 

O ^ 

Q+G{x,-a,,X2) = ^ ^ 

l= — oo 

o 

Q.G{x^+as,X2) = ^ 

/ = —OO 
ID ^ 

PG{xi,X2Tas) = ^ 

l= — oo 

Here Q± are the amplitudes of the two sources located on the xi-axis at 0), whereas P is the amplitude 

of both sources located on the X 2 -axis at (0, Tag); note that these amplitudes are to be found. 

The n-th order coefficient for the total wave incident on the coated inclusion is thus 
Al^) = *" + + ^if(i)(/3ea,)cos(n7r/2). 

As in [18] we set the n-th order coefficient of the Hn^ term equal to zero when r > as and obtain 

4i7o = P Q- P‘^P)['ToHQ^\Pe(^s) P Jo{Pe(^s)]^ 

4ri = {Q^-Q-)[riH[^\peas)PJi{Peas)]. 

4^72 = (Q+ + Q- + 2P) [757^2^^ (/^eUs) + J2(^e<^s)]: 

which we solve for Q± and P for given Finally, we write the displacement outside the coated inclusion 

as 


w{r,0) ~ n;o(xi)-1-(5_ G(r cos^-1-Ug, r sin-j-(5+G(r cos^ — Ug, r sin 

N 

+ P G{r cos r sin 0 — ag) P P G{r cos r sin 0 P ag) P 

n=-N 

where N needs to be chosen sufficiently large to ensure wave amplitude accuracy {N < 8 for all the examples 
considered in this paper). 

In Fig. [Toj^a), we present the membrane wave amplitude for the frequency uj = 2.62 for an inclusion 
with a suitably designed coating chosen to give the flattened variation of |Po^^| in Fig. [^a); this is the 
frequency value at which a small deviation creates large fluctuations in the monopole term. Gomparing 
Fig. [Toj^a) with Fig. |5(b)| (left) or even with Fig. |5(c)| (left), it is clear that the increase in strength of the 
monopole term results in higher scattering. Figures^lop)) and (c) illustrate two attempts, with four and eight 
active sources located at 1.57 units away from the origin, to cloak the coated inclusion. The equations to be 
solved to obtain the amplitudes of eight and twelve sources can be found in the Supplementary Material. 
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Figure 10: Membrane wave scattering: (a) Total displacement field exterior to a coated inclusion with 
parameter values uj = 2.62, pi = 1.5, pc = 0.81, pe = 1.0, pi = 0.1, pc = 0.11, pe = 1.0, Uc = 0.77, = 0.5 
over a square region of side length 20. (b) Total displacement field exterior to a coated inclusion in the 
presence of 4 active sources located on the axes (black dots), all 1.57 units away from the origin. Parameter 
values as in (a), (c) Total displacement field exterior to a coated inclusion in the presence of 8 active sources 
(black dots), all 1.57 units away from the origin; sources off the xi- and X 2 -axes lie on X 2 = ±xi. Parameter 
values as in (a). 


We also remark that employing four active sources is perfectly adequate for smaller frequencies. To 
illustrate this, we present cloaking attempts at frequency values uj = 0.7, 0.8,1.2 denoted by t,t, ft in Fig. 
[Ha), respectively. These frequency values respectively correspond to small, large and intermediate values 
of the monopole term, as shown in the same figure using dotted lines. The contour plots for the membrane 
wave amplitudes when the coated inclusion is surrounded by four active control sources are presented in 
Fig. (a) corresponds to cj = 0.7, (b) to uj = 0.8 and (c) to cj = 1.2. In each of these plots, that on 
the left is with no sources, whereas the middle and the right ones have four sources, the last showing the 
membrane wave amplitude over a wider range of xi^X 2 values. These computations confirm that the values 
of frequencies in the neighbourhood of the points of maximum for (see Fig. correspond to highly 

localised vibration modes, as seen in Fig. |ll(b)[ 

For the higher frequencies of Fig. it is necessary to use far more than four sources to achieve 
good quality cloaking. This is demonstrated in Fig. where cloaking of very good accuracy is achieved 
using twelve sources. This observation is in good accord with the plot given in Fig. of the magnitudes 
of multipole field components. Bearing in mind that the number of active sources is twice the value of the 
order showing on the horizontal axis of Fig. it is evident that good quality cloaking requires ten or more 
active sources. 
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Figure 11: Membrane wave scattering: total displacement field exterior to a coated inclusion, where material 
properties are as in Fig. 10, for (a) uj = 0.7, (b) uj = 0.8, (c) uj = 1.2. Left: uncloaked coated inclusion; 


centre: coated inclusion cloaked using four active control sources (small black dots surrounding the coated 
inclusion); right: same as the central images but over a wider range of Xi, X 2 values. 
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Figure 12: Total displacement field for membrane wave scattering in the presence of 12 active sources (small 
black dots surrounding the inclusion) all 1.57 units away from the origin: (a) exterior to the coated inclusion, 
(b) the entire plate including the inclusion and the coating, (c) as in (b) but a detailed view. The material 


properties and frequency value are the same as in Fig. 10 


Table shows the total outgoing wave coefficients propagating away from the configuration of active 
sources and coated inclusion. For the case where there are no active sources, we see that there is a core set 
of multipole orders in n G [—2, 2] with magnitudes of order unity. When active sources are introduced, it 
is clear that the method reduces to zero the targeted coefficients (see shaded cells). We see an increase in 
magnitude for non-targeted coefhcents which is generally significant for the two orders adjacent to the core. 
Noting for example, the magnitude 0.15 for n = ±6 with eight sources we see why the cloaking in Fig. |lQ(c) 
is still not adequate, but with twelve sources, we see from Table that it is clearly adequate. 



n 

no sources 

4 sources 

8 sources 

12 sources 


-8 

-1.34 X 10-^« + 1.16 X 10-« i 

4.14 X 10“^^ - 0.036 z 

4.14 X 10-16 - 0.036Z 

3.01 X 10-1^ - 0.00026 z 


-7 

-6.25 X 10-'^ - 3.91 X 

0.043 + 2.67 X 10-8 i 

0.043 + 2.67 X 10-« i 

-0.000092 - 5.72 x lO-n^ 


-6 

6.41 X 10“^*^ - 0.000025 z 

-3.68 X lO-'’ + 0.15Z 

-3.68 X 10-6 + 0.15 z 

-9.72 X 10-1^+3.52 x lO-^^z 


-5 

0.00074 + 5.50 X 10“^ i 

0.35 +0.00026 i 

0.023 +0.000017 i 

1.53 X 10-11^ + 3.13 X lO-ii'i 


-4 

-0.00026 +0.016 i 

0.036 - 2.21 i 

-8.45 X lO-ii' - 5.90 X lO-^^i 

-8.38 X 10-1^ + 1.08 X 10-16 i 


-3 

-0.066-0.0044 i 

0.99 +0.065 i 

8.3 X 10-1^ - 4.94 X lO-ii' i 

-1.11 X 10-16 + 1.47 X 10-16 i 


-2 

0.218-0.413Z 

0 

3.04 X 10-16 +3.33 X 10-16 i 

1.25 X 10-16 


-1 

0.305 + 0.896Z 

-1.11 X 10-16 

1.11 X 10-16 - 3.47 X 10-1^ z 

-2.22 X 10-16 + 2.43 x lO-H" i 

0 

-0.963-0.188i 

-5.55 X lO-ii' - 4.44 x lO”!® i 

-5.55 X 10-1^ + 4.44 x 10-16 i 

-1.67 X 10-16 + 4.44 X 10-16 i 


1 

-0.305 - 0.896 

1.11 X 10-16 

-1.11 X 10-16 + 3.47 X lO-ii'i 

2.22 X 10-16 - 2.43 x lO-H'z 


2 

0.218-0.413Z 

0 

3.04 X 10-16 +3.33 X 10-16 i 

1.25 X 10-16 


3 

0.066 +0.0044 z 

-0.99-0.065 z 

-8.3 X lO-ii' + 4.94 X 10-1^ i 

1.11 X 10-16 - 1.47 X 10-16 z 


4 

-0.00026 + 0.016Z 

0.036 - 2.21 i 

-8.45 X lO-ii' - 5.90 X lO-H'i 

-8.38 X 10-1^ + 1.08 X 10-16 i 


5 

-0.00074 - 5.50 X 

-0.35 - 0.00026 z 

-0.023 - 0.000017 z 

-1.53 X 10-1^ - 3.13 X lO-ii'z 


6 

6.41 X 10-^'=’ - 0.000025 i 

-3.68 X 10-6 + 0.15Z 

-3.68 X 10-6 + 0.15 z 

-9.72 X 10-1^ + 3.52 x lO-^^i 


7 

6.25 X 10“"'' + 3.91 X 10“^3 i 

-0.043- 2.67 X 10-8 z 

0.043 + 2.6 X 10-«z 

0.000092 + 5.72 x lO-n z 


8 

-1.34 X 10“^^ + 1.16 X 10-» i 

4.14 X 10-1^ - 0.036 z 

4.14 X 10-1^ - 0.036Z 

3.01 X 10-1^ - 0.00026 z 


Table 1: The coefficients En^ of Hn\Pe'r’) terms, for a configuration with zero, four, eight and twelve control 
sources positioned symmetrically on a circle of radius of 1.57 units from the origin. The frequency of the 
incident wave is uj = 2.62. 


4 Flexural waves for Kirchhoff plates with coated inclusions: res¬ 
onant regimes 


As discussed in Sectional the fields (x), fc = i, c, e, of (2.2) are represented as a sum of solutions of 
the Helmholtz and modified Helmholtz equations, and respectively (see (2.3), (2.4)). In particular, 
when mass densities differ from the matrix to the inclusion, but elastic properties, flexural rigidities and 
Poisson’s ratio, of the inclusion, coating and ambient matrix are the same, the four transmission conditions 
(2.6)-(2.8) on each circular interface reduce to conditions of continuity of the displacement w and its three 
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derivatives with respect to r, i.e. in terms of jump conditions 
gk 

[re(x)] = 0, [^-^re(x)] = 0, with /c = 1, 2, 3, and x being a point on a circular interface. (4.1) 

In this case and at large wavelengths (compared with the inclusion size) the displacement in the matrix is 
dominated by the Helmholtz wave, which satisfies the first two interface conditions (2.6) for the displacement 


and its normal derivative, with the discrepancy in two remaining interface conditions being asymptotically 


small (see Supplementary Material). The mass balance formula (1.1), applied for membrane waves to make 
the monopole term in the scattered field vanish, can then be used as an approximation. However, in the 
intermediate regime when the frequency of the incident wave becomes sufficiently large, the coupling between 
and becomes essential through the interface conditions (|2.7|), (|2.8|). In the intermediate regime. 


the formula (1.1) is no longer valid, just as it is invalid when the elastic properties of the inclusion, coating 


and ambient matrix are different. 

As shown in the previous section for membrane waves, the coating around an inclusion can be used 
to control the monopole term in the expansion of the outgoing wave. Note also that the desired control is 
not necessarily aimed at making the monopole term vanish. Near the resonant frequencies of membranes, 
we have shown that the active cloaking cannot be achieved using the simple mass balance argument of m- 
On the other hand, the use of the coating enables one to shift away from a resonance mode, and hence to 
implement high-precision active cloaking by multiple sources. 

Here, we show that Kirchhoff plates allow a similar treatment of resonant regimes. However, the 
coated inclusion argument requires a significant rethinking, as the mass balance argument, working very well 
for membrane waves, proves to be inadequate for the case of Kirchhoff plates. As expected, the additional 
features of flexural Kirchhoff waves, scattered by an inclusion, are rel ated to the interaction between the 


Helmholtz and the modified Helmholtz waves and defined in (|2.4|). 


(k) 


4.1 Monopole term due to Konenkov 


The scattering of a monochromatic plane flexural wave by a circular obstacle with various boundary 
conditions was intensively studied by Konenkov m- In section 4 of ED, the author gives a formula for 
the scattered field for small values of kR (in original notations k stands for the wave number and R for the 
radius of the inclusion; this corresponds to /3eai in our notations) at a large distance away from the obstacle. 
This formula translates to 


W(sc) ^ iH'^\per 


Pi. 

Pe 

^(1 ~ P'e) + (1 + ^i) 

+ 2 

%)^Ve - (1 + Vi) 

De 

Di 

(1 + Z/^) + ^(1 — Ue) 



{Peaif 


Pedi 


0 . 


With the aim of making the monopole coefficient asymptotically zero to the leading-order as p^ai —> 0, 
we require 


Pei ^e-) Dq) — 

Pe 


^(1 - Ve) + (1 + Vi) 


n 


Vi) 


= 0 . 


(4.2) 


This immediately highlights the difference between membrane and flexural waves in our investigation: all 


material properties of the inclusion and the exterior are present in equation (4.2) for flexural waves, whereas 


for membrane waves the vanishing of to the leading-order merely involves the densities (see formula 


As in the case of an uncoated inclusion in the presence of membrane waves (section 3.2), a strong 


resonant behaviour of the scattering coefficient for a high-contrast inclusion is observed for flexural 
waves. Figure 13^a) shows \E^\Pq)\ for a specific choice of high-contrast material properties pi = 0.05, 
pe = 1.0, Di = 1.5 X 10“^, De = 1.0 and = 0.50. In Fig. 13T) we illustrate the behaviour of \E^\l3e)\ 


for small uj (solid blue rectangle in Fig. p^a)). We see the quadratic behaviour of the monopole term for 
small pe'i this is expected according to Konenkov’s asymptotic formula (see the top right formula on p.l56 
of [21]) in which he states this behaviour is 


^0 


^i^Pii Pei ^ii ^ei ^ii -^e) (/^e^i) ^S Pe^i ^ 0? 
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(in the original notations J(z^, 5, cr, a') {kR)‘^, kR -> 0). 


To clarify, it is important to note that a comparison between Fig. 13 and Fig. 


I requires IF’, 


(e)| 


to 


be plotted as a function of /3e for flexural waves; this is because the spectral parameters are different for 
membrane wave and flexural wave problems, that is jSk = Pk/l^k and Pi = UJ^J pkh/ respectively. 
In other words, |Fq^^| versus uj for membrane waves corresponds to |Fq^^| versus Pe (up to a multiplicative 
constant which is directly comparable to Fig. 


IF, 


(e). 
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Figure 13: Flexural wave scattering by an uncoated inclusion with the parameters pi = 0.05, pe = 1.0, 
Di = 2.5 X 10“^, Fe = 1.0, ai = 0.50: (a) |Fq^^| as a function of Pe (lower horizontal axis) and of uj (upper 
horizontal axis), (b) Detail from (a) showing |Fq^^| (green/dashed) and the empirical fit 0.086 Pi (grey/solid). 


For three different frequencies we give flexural wave spatial distributions in Fig. For the low 

frequency {uj = 0.3) the scattering amplitudes are small. The second frequency considered is uj = 11.15, 
the first non-trivial zero of |Fq^^| where again vibrational amplitudes are small. Going from uj = 11.15 to 
UJ = 11.30, causes dramatic changes in both the spatial distribution around the inclusion and in its amplitude 
variation. 
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Figure 14: Flexural wave scattering by an uncoated inclusion for the parameters as in Fig. 13 (a) Total 
field for uj = 0.3, (b) Total field for uj = 11.15 where vanishes, (c) Total field for uj = 11.30 where 
has a relatively large value. On the right, in each image we present a zoomed-in version of the displacement 
inside and in the close vicinity of the inclusion. The first two frames in the inset (b) represent the same 
field outside the inclusion, but with different colour maps, where the first frame excludes the internal field, 
whereas the second frame includes it. The same comment applies to the first two frames of the inset (c). In 
all contour plots horizontal and vertical axes are the xi- and X 2 -axes, respectively. 


4.2 Coated inclusion: control of the monopole term in the scattered field 

To obtain a uniform asymptotic formula for in the case of a coated inclusion in a thin plate is a 
daunting task. Nonetheless, our aim in this section is to make the monopole coefficient Eq^^ asymptotically 
zero to the leading order, which compared to the previous is somewhat easier. The equation which must to 
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be satisfied in order to make the scattering coefficient zero (to the leadingorder) can be derived as 


[2Z)eai(bi + 2D, - 2D,)^t + - 2D,aihi^t\ «c + [2£»c/3c - 2Aa2b2/3f] 4 

+[4Aa2(b2 + De)4 + a2bl4cPl - 2Aaibi/3,^ + 2A(aibi - a2b2)Pt] 44 = 0, (4.3) 


where 


ai — Di{l 4-I'i) 4-Dc{l — i^c)^ ci2 — Di{l 4-I'i) — Dc{l 4- 

bi = Dc{l ^ l^c) ^ De{l - l^e), b2 = Dc{l - l^c) - De{l - l^e)- 


This fairly unpleasant formula should of course lead us to the result obtained by Konenkov m under 
the assumption of identical material parameter values, for example of the exterior and the coating, hence 
reducing the problem to the case when there is no coating around the inclusion. Assuming Dc = De^Uc = 
Ae, pc = pe and thus the left-hand side of equation (4.3) indeed simplifies to Pe, A, ^e) 


in (4.2). 


It is also important to note that in the very special case of Di = Dc = Ui = Uc = Ve but 
pi pc pei equation (4.3) reduces to equation o> which is exactly what is expected. The mass- 


compensation equation (1.1) is, in general, inadequate to annul the leading-order part of the monopole term 




If one uses the mass-compensation formula unthinkingly, the result may be beneficial or detrimental 
but in either case the result will, in general, still be of second order in /^g. Finally, we draw attention to 


equation (4.3) noting that several parameters associated with the coating may be fine tuned to satisfy the 


equation. This is in stark contrast to the case of membrane waves, where density of the coating was the only 
parameter available for fine-tuning to minimise the monopole term. 




0 . 0.25 1 . 2.25 4 . 6.25 9 . 12.25 16 . 20.25 25 . 



\e4\ 


7.84 9. 10.24 11.56 12.96 14.44 



Figure 15: (a) Flexural wave scattering by (i) an uncoated inclusion with parameter values: pi = 0.05, 
pc = 1.0, Vi = Vc — 0-3, Dc = 1.0, Di = 2.5 x 10“^, h = 1.0, = 0.50 (green/dashed), (ii) a coated inclusion 

with parameter values as in (i) but with material properties of the coating reading: pc = 0.005, Vc = 0-3, 
Dc = 2.5 X 10“^ and Ug = 0.77 (blue/solid), (b) Flexural wave scattering for the parameters of (i) and (ii) in 
(a), but for a range of /3g highlighted by the red box in (a). We see strong variation in \ for an uncoated 
inclusion in this range for ^dg, but it is clear that with an appropriate choice of parameter values, we can 
create a region where remains fairly flat. 
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Figure 16: Flexural wave scattering: (a) Total displacement field exterior to a coated inclusion with parameter 
values as in Fig. 15 and cj = 11.15. (b) Total displacement field of a coated inclusion with parameter values 


as in (a) plotted for the whole plate, including the coating and the inclusion. Note that the displacements 
inside the coating and the inclusion dominate in this case and thus the detail in the exterior is more difficult 
to observe, (c) A zoomed-in version of the displacement inside the inclusion, the coating and in the close 
vicinity of the coating. 





Figure 17: Flexural wave scattering: absolute value of the normalised outgoing wave components 
En^Hn \f3eas) against the multipole order n for uj = 11.15. Parameter values as for (blue/solid) curve 
in Fig. pi = 0.05, pc = 0.005, pe = 1.0, Ui = Uq = Uc = 0-3, Di = 2.5 x 10“^, Dc = 2.5 x 10“^, Dc = 1.0, 
h = 1.0, ai = 0.5, Uc = 0.77, = 1.57. 
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4.3 Active sources: cloaking of flexural waves in resonance regimes 


Here we implement the algorithm of Section 2.2 for active cloaking in regimes, where the frequency is 
sufficiently high, so that the monopole coefficient may vary rapidly as a function of frequency. In this case, 
the coating is introduced to provide a smooth behaviour of in the required frequency range. We also 
comment on the low-frequency regime already discussed above for the case of membrane waves. Explicit 
analytical representations for amplitudes of active sources are included in the Supplementary Material. 


4.3.1 Transition from the low-frequency regime 


As noted in the text above, the scattering produced by the elastic inclusion, in the limit when the 
wavelength of the incident plane wave tends to infinity compared to the size of the inclusion, is very weak. 
In such a regime, the scattering pattern is also dominated by the Helmholtz type waves discussed earlier. 


However, as the frequency of the incident wave continues to increase, the role played by the modified 
Helmholtz waves becomes more important, and we show in this section that a small number of active sources 
is capable of producing a substantial reduction of the scattered field. 


In Fig. |18(a)[ we present the field plot, where uj = 1.50, and hence = 1.2247, the wavelength of the 
incident plane wave is evidently comparable with the diameter of the elastic inclusion. The geometrical and 
elastic parameters of the inclusion and the coating are the same as in Fig. 16 Figure 18(a)| shows a non- 


four active sources 


negligible perturbation of the incident wave by the inclusion. As described in Section [2^ , 
can be introduced to compensate for the action produced by the coated inclusion on the incident field. The 
symmetric configuration of sources, shown in Fig. |18(b)[ is used in the simulation, and the amplitudes of the 
active sources (as outlined in Section 2.2) are derived in an explicit closed form (see Supplementary Material). 
The cloaking with four sources, as in Fig. |18(b)[ is quite successful in compensating for the scattering by 
the coated inclusion for uj = 1.50; beyond this frequency four sources will no longer be adequate. 


In the following section, we introduce larger number of sources and demonstrate how the larger-order 
multipole terms can be suppressed in higher-frequency regimes. 




Figure 18: Flexural wave scattering: (a) Total displacement field exterior to a coated inclusion with material 
properties as in Fig. but at the smaller frequency oi uj = 1.50. (b) Total displacement field for the 
parameter values as in (a) but in the presence of 4 active sources located on the axes (black dots), all 1.57 
units away from the origin. 


4.3.2 Higher-frequency regime 


The multipole representation (2.24) for the total field is used to produce the numerical data discussed 
below. The sources are placed on the circle of radius = 1.57, outside the coated inclusion, whose external 
radius is Uc = 0.77. The simulations are presented for the cases of 4, 8 and 12 sources, together with tabular 
data on the values of high-order multipole coefficients. 
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In Fig. we present three cases corresponding to an incident plane wave scattered by a coated 
inclusion. The parameters of the coating are chosen to be the same as in the simulation depicted in Figure 
16 The case when active sources are absent (Fig. 19 r) is compared with the two cases incorporating four 
and eight active sources, respectively. Although the scattering pattern has been reduced, in this higher- 
frequency regime it is clear from Fig. [^that around twelve sources are required to make the significant 
multipole coefficients vanish. The required combination of twelve sources, with the amplitudes computed 
according to (2.25) is given in Fig. 20, which shows exemplary quality of cloaking, with the scattered field 
being suppressed. 



Figure 19: Flexural wave scattering: (a) Total displacement field exterior to a coated inclusion with parameter 
values as in Fig. (b) Total displacement field in the presence of 4 active sources located on the axes 
(black dots), 1.57 units away from the origin, (c) Total displacement field in the presence of 8 active sources 
located symmetrically 1.57 units away from the origin. 



Figure 20: Flexural wave scattering: (a) Total displacement field exterior to a coated inclusion with parameter 
values as in Fig. (b) Total displacement field in the presence of 12 active sources positioned 1.57 units 
away from the origin, (c) Same as (b) but over a smaller range of the plate detailing localised fields in the 
inclusion. 


It is also noted that Figs. 


20 r and 20 3 use a different colour code scheme, so that the diagram (a) 
plots the field in the exterior of the ambient matrix only, while the diagrams (b) and (c) show the field in 
the inclusion, which also indicates the field localisation. 


The data in Table for multipole coefficients demonstrates that the core region of multipole orders 
runs from n G [—3,3]. Our procedure reduces all the amplitudes of targeted orders to zero, while increasing 
the moduli of adjacent orders. By inspecting the amplitudes of orders ±4 for four sources and ±5 for eight 
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sources, we can see that in both cases we have an insufficient number of sources for successful cloaking. For 
twelve sources, both the entries in Table and the displacement field components shown in Fig. confirm 
good quality cloaking. 



n 

no sources 

4 sources 

8 sources 

12 sources 


-8 

-3.98 X 10“^^ + 2.00 X i 

1.50 X 10-^ - 0.075 z 

1.49 X 10-^ - 0.075 z 

7.66 X 10-y - 0.0038 z 


-7 

-0.000048 - 2.31 X i 

-0.015 - 7.16 X 10-'^ i 

-0.015 - 7.16 X lO-'T 

-0.0022 - 1.05 X 10-'^ i 


-6 

6.70 X 10- ' - 0.00082 z 

0.00027-0.33 z 

0.00027-0.33 z 

-3.43 X 10-1^ + 5.55 x lO-i^ i 


-5 

0.0092 +0.000085 i 

-0.055 - 0.00051 z 

0.24 +0.0022 z 

-1.39 X 10-1^ - 8.02 X lO-ii'i 


-4 

-0.0039 +0.0629 z 

0.072 - 1.15 z 

-1.96 X 10-1^ + 1.11 X 10-16 i 

-4.84 X 10-1® - 1.11 X 10-16 z 


-3 

-0.20 - 0.043 z 

-0.26 -0.055 i 

-6.94 X lO-ii' + 1.08 X 10-16 i 

1.39 X lO-ii' - 1.71 X 10-1® z 


-2 

0.066 - 0.25 z 

-9.26 X 10-1® + 8.33 x lO-i^^ i 

-4.87 X 10-1^ + 1.53 X 10-16 i 

4.49 X 10-1® + 1.25 X 10-16 i 


-1 

-0.071 + 0.0051 i 

-1.39 X lO-ii' + 4.50 X 10-1® i 

-4.16 X 10-1^ + 8.71 X lO-ii' i 

-4.16 X 10-1^ + 6.31 X 10-1^ i 

0 

-0.98 - 0.13z 

1.58 X 10-16 - 5.55 X 10-1^ i 

-7.88 X 10-1^ 

-9.27 X lO-ii' - 2.78 x 10-i^z 


1 

0.071 - 0.0051 z 

1.39 X 10-1^-4.50 X 10-1® 

4.16 X lO-ii' - 8.71 X 10-i^i 

4.16 X lO-ii' - 6.31 X lO-ii'i 


2 

0.066 - 0.25 z 

-9.26 X 10-1® + 8.33 x lO-H' 

-4.87 X lO-ii' + 1.53 X 10-16 i 

4.49 X 10-1® + 1.25 X 10-16 i 


3 

0.20 +0.043 z 

0.26 +0.055 z 

6.94 X lO-ii' - 1.08 X 10-16 i 

-1.39 X lO-ii' + 1.71 X 10-1® i 


4 

-0.0039 +0.0629 i 

0.072 - 1.15z 

-1.96 X lO-ii' + 1.11 X 10-16 i 

-4.84 X 10-1® - 1.11 X 10-16 i 


5 

-0.0092 - 0.000085 z 

0.055 +0.00051 z 

-0.24 - 0.0022 z 

1.39 X lO-ii" + 8.02 X lO-ii' i 


6 

6.70 X 10- ' - 0.00082 i 

0.00027-0.33 i 

0.00027-0.33 i 

-3.43 X 10-1^ + 5.55 x i 


7 

0.000048 + 2.31 X i 

0.015 + 7.16 X lO-'T 

0.015 + 7.16 X lO-'T 

0.0022 + 1.05 X 10-'^ i 


8 

-3.98 X 10-^^ + 2.00 X lO-'^ i 

1.50 X 10-1 - 0.075 z 

1.49 X 10-1 - 0.075 z 

7.66 X lO-"" - 0.0038 z 


Table 2: The coefficients En^ of Hn\Pe'f') terms for a configuration with zero, four, eight and twelve control 
sources positioned symmetrically on a circle of radius of 1.57 units from the origin. The frequency of the 
incident plane wave is uj = 11.15. 


5 Concluding remarks 

This paper has presented a novel idea of active cloaking for coated inclusions subjected to incident 
waves of frequencies ranging from zero up to higher frequency regimes, where the monopole term in the 
outgoing scattered field can be a rapidly varying function of frequency. The technique involves choosing the 
coating to diminish the frequency sensitivity of the scattering problem, and then applying active control 
sources. 

A comparison between the membrane waves and flexural waves in Kirchhoff plates has revealed in¬ 
teresting connections for low frequencies. In particular, we have identified the cases when the flexural wave 
scattering is driven by the Helmholtz waves. On the other hand, we have also identified important differences 
where the evanescent waves which are present in the expressions for the solutions to the fourth-order problem 
contribute significantly. 

The utility of the proposed algorithm is high, as the amplitudes of active sources are given in the 
explicit analytical form, and the accuracy of the cloaking of coated inclusions in resonant regimes has been 
well demonstrated. 
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